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Abstract 

At the very heart of the successful phenomenological model of magnetic hysteresis 
there is the so called Preisach distribution. In the existing literature it is implicitly as- 
sumed, that this distribution has a mirror symmetry. We show, by simple and convincing 
example, that this common assumption is plainly wrong. Dropping it, we gain the ability 
to model not only the usual hysteresis loops (major and minor) more accurately than ever 
before, but also those displaying the exchange bias effect, what is impossible within the 
framework of the symmetrical Preisach model. It is hoped, that our observation paves 
the way towards the unified description of all hysteretic systems, including, but not nec- 
essarily limited to, superconductors, (multi)laycrcd structures, nanocrystalline materials, 
patterned media, and - perhaps - the other non-magnetic hysteretic phenomena. 

Introduction The major hysteresis loop, observed in the sizable samples of homogeneous 
ferromagnetic materials, exhibits well known symmetry: 

Mi b (ff) = -Mub(-JJ), (1) 

where Mib (M u b) denotes the lower (upper) branch of the sample's magnetization, M vs. 
exciting field H . The hysteresis curve, not necessarily the major loop, but also the response 
to the arbitrary sequence of exciting fields as well, can be described, or modelled, in many 
ways. One of them, the Classical Preisach Model (CPM), was first proposed by Ferenc 
Preisach [1] and then subsequently developed, generalized and tested by many researchers. 
In this model, the change of the sample's magnetization is expressed by the double integral: 

AM = M(„,)-M<*,= 2 M. jj ,(„„*). <2> 

H i< H l< H { 

for the monotonously increasing field ('i' - the initial state, 'f - final) and by 

AM = M(H { ) - M(Hi) = -2M S / J dff T dH l g(H h Hi), (3) 

h t >s 1 

H f< H X < H i 

if the field is monotonously decreasing. M s is the saturation magnetization, while the distri- 
bution g(Hf,Hi) > 0, supported over the domain > Hi is called the Preisach density. 
With the additional condition: 

dff T dffj. q(H u H 1 ) = 1, (4) 
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the Preisach density is commonly regarded as a probability distribution. It describes the 
probability of encountering the so called hysteron, i.e. the hypothetical object characterized 
by elementary, rectangularly shaped, hysteresis loop with up- and down-switching fields equal 
to Hi and Hy, respectively. 

It is obvious, that the symmetry presented in Eq. (1) will be preserved, if the Preisach 
distribution, g, is also symmetric, namely when 



g{H t ,H l ) = Q{-H u -Hi 



(5) 



holds for any pair (-Hj- , H±) of its arguments, Hf > Hy. The mirror symmetry of g (Hf, H\) 
with respect to the line = —Hf is therefore commonly assumed, before any reconstruction 
of the Preisach distribution is ever attempted. In this paper we are going to show, that 
symmetric hysteresis loops can be generated by Preisach distributions, which themselves 
have no mirror symmetry. In other words: while the mirror symmetry of g indeed implies 




Figure 1: An example of the non-symmetric Preisach map, which generates symmetric 
major hysteresis loop. Initially the only components of the map are two symmetrically 
located pairs of positive, delta-shaped peaks, with identical amplitudes, labeled as (P, Q) 
and (P',Q') respectively 



the experimentally observed symmetry properties of hysteresis loops, then the inverse need 
not to be true. The idea that the Preisach distribution has to be symmetric is present in the 
literature since almost 50 years [2]. It seems so obvious that nobody ever cared to prove it 
formally. For example, Mayergoyz (1986) wrote in [3] 

... it can be easily proved that u, (a, (3) = n (—(3, —a) 

however he didn't present the proof (Mayergoyz's n is our g). 



Hints for possible asymmetry Let us try to prove it now by reduction ad absurdum. 
Assume that we have a nonsymmetric distribution and let's see whether it necessarily has 
to generate the non-symmetric hysteresis loop. The distribution is non-symmetric, if there 
exists at least one point (H^ , H*^ in its support, such that the condition g (//| , J/j) = 
g(—H*,—H*^ is violated. But if this indeed is the case, then the integrals (2) and (3) 
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remain unchanged, since a single point is a zero-measure set on a Preisach plane! As we can 
see, it is impossible to prove the well known claim. One could argue that our finding has no 
physical meaning, since such kind of asymmetry cannot be detected in any real experiment 
and, consequently, it should be regarded only as a mathematical toy. Eventually we might 
modify the original claim to read: the Preisach distribution is symmetric almost everywhere 
on a Preisach plane. However, we are not going to stop at this point and shall present 
a well-founded, physical counterexample. 
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Figure 2: Stacked hysteresis loops generated by non-symmetric Preisach distribution pre- 
sented in Fig. 1 for various amplitudes of the exciting field [-H, +H]. For small amplitudes 
there is a null response (hence the loop is symmetric), in the intermediate regime the loop 
becomes asymmetric and exhibits the so called exchange bias effect, while the major loop 
(lowest) is symmetric again 

Derivation of the main result Our counterexample uses four Dirac's delta functions as 
the sole components of a Preisach distribution. This is not completely crazy idea, since the 
delta peaks are the only mathematical entities able to reproduce well known Barkhausen 
jumps. Those jumps appear sometimes very stable against repeated magnetization reversals. 
The same is true for the domain structures [4, 5]. The other advantage is that the single-point 
support of a delta function is no longer a zero-mesure set for integrals (2) and (3). We are 
not original introducing delta functions into Preisach model (see Pescetti [6] 1989), so our 
derivation cannot be viewed only as a mathematical curiosity. 

Consider Fig. 1. Initially the Preisach distribution consists of two delta peaks of equal 
amplitudes, located at points labeled as P and Q, and their symmetrical images P' and Q' , 
respectively. Such a distribution gives rise, of course, to symmetrical hysteresis loops. But 
now we destroy the initial symmetry by moving points Q — ► Q" and P' — * P", as indicated 
by arrows. Having done so, we can observe (Fig. 2) that the major hysteresis loop, as well 
as some minor loops, remain symmetric, although their overall shape has changed. There 
are also some minor loops, which are asymmetric now, imitating the so called exchange bias 
effect. 

It is also possible to rearrange the four initial delta peaks differently, in such a way as to 
preserve the symmetry of all minor loops with turning points of equal magnitudes. This can 
be achieved by rotating the pairs (P, Q) and (P 1 , Q') around points P and Q' respectively by 
the same angle ip. Let initialy P — (Ho, —H2), Q = (Ho, —Hi), with d = \Hi — Hi\. After 
rotation, the jumps on the lower branch of the hysteresis loop ccur at fields: Ho, Ho + d sin <p, 
Hi and Hi + dcosip (in increasing order). Upper branch jumps (in decreasing order) are 
located at: — Ho, —Ho — dsin ip, —H2 + dcos <p and —H2. One can easily see, that in all cases 
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i?! J ™ p = —H 1 v ^ np , what proves that all the relevant symmetry properties of the hysteresis 
loops are indeed preserved. 

Comments Our construction strongly suggests that the true (asymmetric) Preisach dis- 
tribution may be a conformal image of some 'ideal', i.e. symmetric map. It is worth to 
note, however, that conformal mapping of the Preisach triangle onto itself, if applicable at 
all, cannot be arbitrary. Ordinary Euclidean transformations, i.e. translations and rotations, 
are useless (shifted distribution fails to describe exchange biased loops, see [7]). It remains 
unclear, whether the asymmetry of the Preisach distribution alone would be able to describe 
correctly the recently reported [8] hysteresis loops exhibiting negative remanence. 

Conclusions We have shown that, contrary to the common belief, the Preisach distribution 
need not to be symmetric. This seems to remove the apparent disagreement between the 
distributions recovered from experimental data on one of the classical ways [9, 10, 11] (which 
implicitly assume the existence of symmetry questioned here, sometimes even including the 
postulated particular shape of the map to be reconstructed) with those obtained by the 
recently introduced FORC (First Order Reversal Curves) [12, 13] diagrams technique, which 
reveals the asymmetry quite often; see [14, 15, 16] and especially impressive figures presented 
by Robb, Novotny and Rikvold [17]. 

The construction of FORC diagrams implicitly suppresses the reversible part of the Preisach 
map, i.e. the one located on the line = H±. In our opinion, it is the combination of both 
techniques, what should constitute the complete tool for physical, as opposed to purely phe- 
nomenological, characterization of magnetic systems and their internal interactions. 
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